
Available online at www.sciencedirect.com 



r* • r*\* ■ Nuclear Physics B 

V ScienceDirect proceedings 



Supplements 

ELSEVIER Nuclear Physics B (Proc. Suppl.) 00 (201 1) l-[5] ^^^^^^^^= 



Measurement of the two-photon exchange contribution 
in elastic ep scattering at VEPP-3 

A.V. Gramolin a '*, J. Arrington b , L.M. Barkov a , V.F. Dmitriev a ' c , V.V. Gauzshtein d , R.A. Golovin a , RJ. Holt b , 
V.V. Kaminsky a , B.A. Lazarenko a , S.I. Mishnev a , N.Yu. Muchnoi a ' c , V.V. Neufeld a , D.M. Nikolenko a , LA. Rachek a , 
R.Sh. Sadykov a , Yu.V Shestakov ac , V.N. Stibunov d , D.K. Toporkov ac , H. de Vries 6 , S.A. Zevakov a , V.N. Zhilich a 

a Budker Institute of Nuclear Physics of SB RAS, Novosibirsk, Russia 
b Argonne National Laboratory, Argonne, USA 
c Novosibirsk State University, Novosibirsk, Russia 
d Nuclear Physics Institute of Tomsk Poly technic University, Tomsk, Russia 
e NIKHEF, Amsterdam, The Netherlands 



Abstract 

We report on the status of the Novosibirsk experiment on a precision measurement of the ratio R of the elastic e + p 
and e~p scattering cross sections. Such measurements determine the two-photon exchange effect in elastic electron- 
proton scattering. The experiment is conducted at the VEPP-3 storage ring using a hydrogen internal gas target. 
The ratio R is measured with a beam energy of 1.6 GeV (electron/positron scattering angles are 6 = 55-^75° and 
6 = 15-^25°) and 1 GeV (0 = 65-^105°). We briefly describe the experimental method, paying special attention to the 
radiative corrections. Some preliminary results are presented. 
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1. Introduction 

The elastic electromagnetic form factors of the pro- 
ton are important characteristics that describe its inter- 
nal structure. More precisely, the electric Ge(Q 2 ) and 
magnetic G M (Q 2 ) form factors describe the spatial dis- 
tributions of the proton's charge and magnetization as 
functions of four-momentum transfer squared Q 2 (see 
Refs. dEl for reviews). Until recently these form fac- 
tors were determined only by an analysis of differential 
cross section of the unpolarized elastic ep scattering. 

Let us recall that the differential cross section dcr/dQ 
for elastic ep scattering in the one-photon exchange ap- 
proximation (see the first diagram in Fig. [TJ is given by 
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the Rosenbluth formula (3) : 

do~ t r 2 2 s r 2 (n 2 \\ d^Moa 

where r = g 2 /(4M 2 ), s = [1 + 2(1 + r)tan 2 (<9/2)] _1 is 
the virtual photon polarization, is the electron scatter- 
ing angle in the laboratory frame and dcr Mott /dQ is the 
Mott differential cross section. It follows from this for- 
mula that the proton form factors can be determined by 
measuring the differential cross section at fixed momen- 
tum transfer g 2 , but with different electron scattering 
angles and incident beam energies. 

In the mid-nineties, it became possible to use another 
method to study the proton form factors, the polariza- 
tion transfer method, which was originally proposed 
back in 1974 [4, 5 ]. In this method a polarization of the 
recoil proton in the process of elastic scattering of lon- 
gitudinally polarized electrons on an unpolarized hydro- 
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Figure 1: Feynman diagrams of the two lowest orders in a for the 
elastic electron-proton scattering. The first diagram corresponds to 
the one-photon exchange approximation. The box and crossed-box 
diagrams describe the two-photon exchange (the blob represents pos- 
sible intermediate states of the proton). 



gen target is measured. In such a case the ratio Ge/Gm 
is directly proportional to the ratio of transverse and lon- 
gitudinal polarizations of recoil protons. 

Measurements performed at TJNAF with the polar- 
ization transfer method |[6lfT0ll gave the results that are 
in contradiction with the results obtained by the Rosen- 
bluth method. This is clearly seen in Fig. [2| where the 
polarized data of TJNAF as well as several results of 
the Rosenbluth measurements are shown. The latter in- 
clude the global analysis by Walker [11] and new accu- 
rate Rosenbluth measurement lfT2l . 

The most probable cause of this discrepancy is the 
failure of the one-photon exchange approximation to 
precisely describe the results of unpolarized measure- 
ments. However, application of the two-photon ex- 
change corrections (see the diagrams in Fig. [I]) runs into 
difficulties: on the one hand, there are no reliable calcu- 
lations and, on the other hand, there are no sufficiently 
accurate experimental data (for review see Ref. fT3ll ). 

The two-photon exchange effect can be measured ex- 
perimentally by precisely comparing the e~p and e + p 
scattering |[T3l . It is possible because the interference 
between the one-photon exchange amplitude and the 
two-photon exchange amplitude occurs with opposite 
sign for electrons and positrons. In 2004, it was pro- 
posed lfT4l to perform such measurement at the VEPP-3 
storage ring H31 in Novosibirsk. Here we report on 
the status of this experiment. It should be mentioned 
that there are also two similar experiments underway at 
DESY (OLYMPUS collaboration OH) and at TJNAF 
(CLAS collaboration Pi ). 



2. Description of the experiment 

The experiment was divided into two runs with dif- 
ferent kinematic parameters (for details see Table [T}. 
The first run (with a beam energy of 1.6 GeV) was con- 
ducted in late 2009. Analysis of data obtained in this 
run is continuing. The preliminary results are presented 
here (also see Ref. fT8ll ). The second run (with a beam 




Figure 2: Comparison of the polarized (oTflOl and unpolarized [TT] 
[T2l data for the ratio juGe/Gm (where \x is the proton magnetic mo- 
ment). The solid lines are the quadratic fits for these two groups of 
data. 



energy of 1 GeV) is underway and will continue, pre- 
sumably, until March 2012. 

The ratio R = cr(e + p) / cr(e~ p) of the elastic positron- 
proton and electron-proton scattering cross sections is 
measured. The scattered electron/positron and the recoil 
proton are detected in coincidence. In the first run, the 
particles were detected in three angular ranges: large an- 
gles (LA), medium angles (MA) and small angles (SA). 
Scattering of e~ /e + at small angles was used only for 
luminosity monitoring (this is possible because the two- 
photon exchange effect is negligible for this kinematic 
range). In the second run, the particles are detected in 
two angular ranges: large angles (LA) and medium an- 
gles (MA). The latter range is used for luminosity mon- 
itoring. Kinematic parameters for both runs are shown 
in Tableffl 

The hydrogen internal gas target is used, similar to 
that previously used in the experiments at VEPP-3 lfT9l . 
The target is an open-ended storage cell, which has an 





Run I 


Run II 


Parameter 


LA 


MA 


SA 


LA 


MA 


^beam, GeV 


1.6 


1.0 


6, degree 


55/ 


15/ 


8/ 


65/ 


15/ 


(min/max) 


75 


25 


15 


105 


25 


G 2 , GeV 2 


1.26/ 


0.16/ 


0.05/ 


0.71/ 


0.07/ 


(min/max) 


1.68 


0.41 


0.16 


1.08 


0.17 


s 


0.58/ 


0.97/ 


0.99/ 


0.51/ 


0.97/ 


(max/min) 


0.37 


0.90 


0.97 


0.18 


0.91 



Table 1 : Kinematic parameters of the experiment. 
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Figure 3: Schematic side view of the detector used in the first run. 



elliptical cross-section 13 x 24 mm 2 and a length of 
400 mm. High-purity hydrogen gas is directed into the 
center of the storage cell. The thickness of the target 
is about 10 15 at./cm 2 . Hydrogen is pumped from the 
VEPP-3 vacuum chamber using cryogenic pumps. 

The detector (see Fig. [3]) consists of two identical 
parts (top and bottom) arranged symmetrically with re- 
spect to the median plane of the storage ring. Azimuthal 
angle acceptance of the detector is Acp = 2 • 60° = 120°. 
In the first run, the particles were detected in three angu- 
lar ranges (LA, MA and SA). In the LA and MA ranges 
the detector has a tracking system (drift and multi-wire 
proportional chambers), calorimeters (Csl and Nal crys- 
tals) and plastic scintillators used for event triggering 
and for time-of-flight measurements. As for the SA 
range, here we used the sandwich scintillators. 

In the second run, we have a slightly different config- 
uration of the detector, but it is similar to that described 
above. An important feature is that here we have the 
possibility to reconstruct the full kinematics for some 
of the events with emission of hard real photons. This 
will allow us to test our calculations for the first-order 
bremsstrahlung (see the next section). 

We use the following cuts (in various combinations 
for different ranges) for the selection of elastic scatter- 
ing events: correlation between the polar angles of e~/e + 
and p\ the same for the azimuthal angles; correlation be- 
tween the lepton scattering angle and the proton energy; 
correlation between the lepton scattering angle and the 
lepton energy; AE-E analysis for the proton identifica- 
tion; and time-of-flight analysis for the identification of 
low-energy protons. The main sources of background 
events are pion production processes and machine back- 
grounds. 



To suppress systematic errors it is important to main- 
tain the same experimental conditions during data tak- 
ing with electron and positron beams. Therefore, data 
collection with e~ and e + beams is alternated regularly. 
This allows us to suppress effects of slow drift in time of 
the target thickness, detection efficiency and some other 
parameters. 

A typical working cycle of the experiment consists 
of the following stages: storing of e~/e + beam up to 
40 -r 55 mA; ramping beam energy (from the injection 
energy of 0.35 GeV to 1.6 or 1 GeV); switching on the 
internal gas target and then data taking (as long as the 
beam current drops to 10 -r 13 mA); and returning to 
the injection energy. The full cycle with two beams 
(electron and positron) takes 60 -r 80 minutes, from 
which the data taking occupies approximately 40 min- 
utes. The average luminosity of the experiment is about 
5 • 10 31 s" 1 • cm -2 . The beam current integral collected 
in the first run is 54 kC. In the second run, we plan to 
collect the beam current integral of 60 -f 90 kC. 

The main possible sources of systematic errors of the 
measurement are as follows: different beam positions 
for electrons and positrons; unequal beam energy for 
electrons and positrons; time instability of detection ef- 
ficiency; drift of the target thickness in time; and uncer- 
tainty related to the radiative corrections. For the first 
run, the total systematic uncertainty is now estimated 
as 3 • 10" 3 (for both LA and MA ranges). 

There are three sources of information on the beam 
position: VEPP-3 beam position monitors (electrostatic 
pick-ups); movable beam scrapers; and reconstruction 
of the event vertex position using the tracking system 
of the detector. As for the beam energy, it is measured 
regularly and with good accuracy (about 50 keV) dur- 
ing the data taking by applying the Compton backscat- 
tering technique (20l ETJ. The essence of this tech- 
nique is as follows. Photons from a CO2 laser are scat- 
tered in a head-on collision with the stored beam. The 
beam energy can be determined from the spectrum of 
the backscattered photons that are detected by a high- 
purity germanium detector. 

3. Radiative corrections 

Not only the two-photon exchange, but also other 
higher-order QED processes affect the elastic scatter- 
ing cross section. These processes include internal 
bremsstrahlung of real photons (see the diagrams in 
Fig. [4]) and the processes associated with virtual photons 
(vacuum polarization and vertex corrections, see Fig. [5}. 
Note that the external bremsstrahlung is neglected due 
to the small thickness of the target. 
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Figure 4: Feynman diagrams representing the first-order internal 
bremsstrahlung. 



Figure 5: Feynman diagrams representing virtual photon corrections. 



An important fact is that the internal bremsstrahlung 
affects the measured ratio R (because the interfer- 
ence between the lepton bremsstrahlung and the proton 
bremsstrahlung occurs with opposite sign for electrons 
and positrons). The bremsstrahlung process can easily 
lead to a difference of a few percent between the cross 
sections of e~p and e + p elastic scattering. It is there- 
fore very important to carefully consider the radiative 
corrections due to bremsstrahlung. 

In addition, the radiative corrections due to brems- 
strahlung depend on the detector geometry, its spatial 
and energy resolutions and the kinematic cuts used in 
event selection. Accounting for all these factors is only 
possible with careful simulation of the detector. For this 
purpose, a new Monte Carlo event generator l22l . called 
ESEPP (Elastic Scattering of Electrons and Positrons 
by Protons), has been developed. The generated events 
(e ± p — > e ± p and e ± p — > e ± py) are used for Geant4 
simulation of the detector response. 

We apply the exact calculation l23l for the first-order 
bremsstrahlung. The main features of this calculation 
are as follows. First, we do not use the peaking approx- 
imation or even the soft photon approximation. Sec- 
ond, we consider both lepton bremsstrahlung and pro- 
ton bremsstrahlung. Third, in the proton line we use the 
full electromagnetic vertex operator (see Ref. [ 24 1) 
for the real photon emission. In addition, we apply the 
calculation l25l for the first-order bremsstrahlung with 
the delta-isobar A(1232) excitation. We use this calcula- 
tion as an estimate of the contribution of the first-order 
bremsstrahlung with the proton excitations in the inter- 
mediate state. A detailed description of these calcula- 
tions and the ESEPP event generator will be published 
elsewhere. 

As for the other processes (vacuum polarization and 
vertex corrections), we use the standard approach de- 
scribed in Refs. l24l l26l . Multiple soft photon emis- 
sion | '24' | is not taken into account, as the estimates show 
that it can be neglected under our experimental condi- 
tions. As usual, only the infrared divergent parts of the 
two-photon exchange amplitude and the proton vertex 
correction are taken into account. These terms exactly 
cancel with the corresponding infrared divergent terms 
from the bremsstrahlung process. 



4. Preliminary results 

Preliminary results of the first run are shown in Fig. [6] 
in comparison with all available world data on the ra- 
tio R (in the range Q 2 < 2 GeV 2 ). It can be seen 
that our results have better accuracy than the data avail- 
able. For the LA range (e = 0.5, Q 2 = 1.43 GeV 2 ) 
the result is R = 1.016 ± 0.011 ± 0.003, and for the 
MA range (e = 0.95, Q 2 = 0.23 GeV 2 ) the result is 
R = 0.9976 ± 0.0009 ± 0.003. Here, the first error is 
statistical and the second is systematic. These results 
are preliminary because some minor corrections have 
not yet been made (for example, corrections related to 
the variation of beam energy and position). Note that 
before applying the radiative corrections, these values 
were R = 1.053 for the LA range and R = 1.0044 for 
the MA range. 

The curves in Fig. [6] show the theoretical predictions 
for the ratio R due to the two-photon exchange [27 1 . 
These data were kindly provided by Prof. Blunden [28 ]. 
The solid and dashed lines correspond to the kinematics 
of the first and the second runs, respectively. If these 
predictions are correct, then the expected size of the 
two-photon exchange effect in the second run should be 
approximately the same as for the LA range in the first 
run. At the same time, the elastic e ± p scattering cross 
sections are larger for the kinematics of the second run. 

The second run of the experiment is underway, and 
we are not ready to present results yet. We plan to get 
results in two or three experimental points in the range 
s = 0.2 -r 0.5 with total errors of about 6 • 10" 3 . 

5. Conclusion 

We have described the status of the Novosibirsk ex- 
periment on the precise comparison of the e~p and e + p 
elastic scattering cross sections. The results of our 
measurement are more accurate than the data currently 
available. The experimental data give the size of the 
two-photon exchange effect in elastic electron-proton 
scattering. This information is important to explain the 
recently observed discrepancy between new polarized 
measurements and old data on the proton electromag- 
netic form factors. 
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Figure 6: Comparison of the preliminary results for the ratio R (radia- 
tively corrected) with the previous measurements (for Q 2 < 2 GeV 2 ): 
(29); □ (30), first experiment; A (30), second experiment; T f3Tl ; 
X (32l : ■ f33l: A (34) : O f35l : • this experiment, run I. The curves 
show the theoretical predictions [27 28 1 for the ratio R due to the two- 
photon exchange. The solid (dashed) line corresponds to the kinemat- 
ics of run I (run II). 

Special attention is given to the radiative corrections, 
in particular, to the first-order bremsstrahlung. We have 
developed the event generator, which is based on the 
accurate calculation of the first-order bremsstrahlung 
beyond the soft-photon approximation. The generated 
events are used for Geant4 simulation of the detector 
response. This allows us to take into account the geome- 
try of the detector, its coordinate and energy resolutions 
and kinematic cuts used in event selection. 

We have presented the preliminary results of the first 
run. The second run of the experiment is underway. It 
is too early to draw final conclusions based on the re- 
sults presented, since the data analysis is still ongoing. 
It can be noted, however, that the preliminary results are 
consistent with the theoretical predictions l27l[28l . 
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